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D ■ We present an expression for the generating function of correlation functions 

in the sine-Gordon integrable field theory on a cylinder, with compact space. This 
is derived from the Destri-De Vega integrable lattice regularization of the theory, 
formulated as an inhomogeneous Heisenberg XXZ spin chain, and from more recent 
qq . advances in the computations of spin form factors in the thermodynamic limit. 
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in 

1 Introduction 

o 

Integrable massive field theories in two dimensions are studied since long time, serving 
as exactly solvable prototypes of strongly interacting theories. Applications to problems 
in condensed matter and in statistical physics have been found over the years (see, 
e.g., [U E]) and range from off-critical statistical field theories to low-energy effective 
descriptions of conduction and magnetism in one dimension. 

The sine-Gordon (SG) field theory is a paradigmatic example of integrable field theory 
because of both its simplicity and the richness of its features. In this paper, we are 
interested in computing SG correlation functions in Euclidean spacetime on a finite 
geometry, i.e. on a cylinder, on which the compactified direction will be regarded as 
space and the infinite one as time. 

On the plane, the problem of computing correlation functions can be formally solved 
by applying the standard form factor approach [3HU|5]. Conversely, for finite size, while 
studies of form factors have been carried on in semiclassical approximation in [6j [7], 
the lack of a simple expression for the vacuum of the theory and its excitations still 
constitutes an obstacle for the method. 

A way of tackling the problem may be that of interchanging the labels of space 
and time and consider the problem in which space is infinite, while the imaginary time is 
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periodic, which amounts to considering the system at finite temperature. With this setup, 
it is possible to use the "infrared" form factors, but a suitable regularization is nonetheless 
necessary. One-point functions can be efficiently computed by series [8j [9] ; however, the 
formalism for two-point functions, perturbative in the exponential corrections in the 
volume, presently relies on a double series expansion containing regularized infinite- 
volume form factors [ID] , which results in cumbersome expressions for the model under 
consideration. 

Going back to the original perspective, it is known [TT] that corrections to the particle 
masses which constitute the spectrum in the infrared limit are exponentially suppressed 
in the size of the system. This fact allowed the proposal [12] of a finite- volume formalism, 
based on the infinite-volume form factors, that is correct up to terms which are expo- 
nentially suppressed in the volume. The method has been applied in [T31 UJ] to soliton 
and breather form factors. An alternative regularization has been proposed in [15J and 
applied to the computation of the dynamical structure factor in prototypical integrable 
field theories, although it is presently unclear how to extend the scheme to more general 
form factors. 

The goal of this paper is that of deriving an exact expression for the generating func- 
tion of correlation functions in finite volume, that includes exponential corrections in 
the size. The framework will be that of lattice integrable regularizations of the SG field 
theory and in particular of the one proposed in [16] . The expression derived will be writ- 
ten in terms of the counting function which solve the Destri-De Vega nonlinear integral 
equation (NLIE) as an expansion over the basis of the exact finite-volume eigenstates 
of the theory. We will moreover take advantage of the algebraic Bethe Ansatz results 
relative to the one-dimensional (inhomogeneous) Heisenberg magnet and in particular of 
the computation of the matrix elements of the magnetization operator in the limit [17] 
in which the number of sites goes to infinity. 

By means of algebraic Bethe ansatz, a determinant representation for the generat- 
ing function has been the goal of [IB]. We believe that the subsequent advances on 
quantum spin chains (such as [TH]) allow more explicit results. In connection with the 
inhomogeneous XXZ spin chain, the one-point functions of primary fields and their de- 
scendants in the sine-Gordon model have been analyzed in the framework of the recently 
explored fermionic structure of the model in [20] , while the computation of form factors 
by separation-of-variables has been tackled in |2TJ recently. 

Section [2] reviews some basic facts about the field theory and its lattice regulariza- 
tion and allows us to fix the notation. The generating function of connected correlation 
functions is defined in Section [3] from the quantum field theory and from the spin chain 
perspectives and its expression is given in Section HI Manipulation of the lattice quan- 
tities is performed in the subsequent Section then the scaling limit is performed in 
Section 
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2 Sine-Gordon and the inhomogeneous XXZ spin 
chain 



The action of the sine-Gordon field theory on a cylinder of radius L is: 

i«W0-^cos(/30) 



oo rL 

dr / dx 



-oo JO 



where \x and /3 are real parameters. In infinite volume, the fundamental excitations are 
known to be the soliton, with mass m = m(fi 2 ) and unit topological charge, and the 
antisoliton, with equal mass and opposite charge. A soliton and an antisoliton can bind 
together and form a breather, and breather scattering can lead to the production of a 
higher mass breathers as well. To be specific, bound states are labeled by integers and 
their number depends on the value of (3. By defining the parameter p = 8 J*_p2 , it is 
possible to distinguish two regimes: a repulsive one, in which only the soliton and the 
antisoliton are present in the spectrum, and an attractive one, in which a number \l/p\ 
of bound states are allowed. 

Correlation functions of the fields on the plane can be in principle generated by the 
knowledge of the form factors of the exponential of <fr, which have been computed in 
[221 IS] , while field form factors have been studied in [U [231 [21] • 

The lattice construction of [16], based on a vertex model, can be formulated by 
means of an inhomogeneous XXZ spin chain, as underlined by the authors. Essentially, 
the fermionic variables defined on the bonds among vertexes, naturally identify spin 
variables via Jordan- Wigner transformation. The alternating sign rapidities which are 
associated to the bonds constitute the inhomogeneities of the XXZ chain, while the 
interaction among fermions of the original work is related to anisotropy parameter of 
the easy-axis interaction and to the coupling in the sine-Gordon Lagrangian. It has 
been shown in [25| [26[ [271 [28] that this formalism is able to reproduce many features 
of the theory in the limit L — > oo (energy spectrum, scattering phases) as well as those 
expected from conformal field theory (operator scaling dimensions) when L — > 0. 

Consider [21] a lattice Heisenberg spin 1/2 chain with 2N sites, characterized by local 
R-matrices of the form: 

1-c . 1 + c , „ . b 

2 

sinhA sinhry 
sinh (ry — A) sinh (27 — A) 

in which A is a complex " Bethe ansatz" rapidity and 7 is the anisotropy parameter, since 
in the homogeneous case is related to the third-component interaction of the spin degrees 
of freedom. It is related to the field theory parameters above via p — - — 1. The indexes 
and n refer to the auxiliary C 2 local Hilbert space and to the one relative to the n— th 
site. 

The monodromy matrix for an inhomogeneous generalization of the XXZ chain can 
be built by means of an auxiliary C 2 space as follows: 

T(A) = e^W-VRoM* - A + *7/2)iViv-i(A + A + z 7 /2) . . . i? ,i(A + A + z 7 /2) (3) 
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in which twisted boundary condition with twist parameter u have been introduced. 

Eigenstates of the transfer can be written through the application of generalized 
creation operators B(X) on a reference state, chosen here to be the one in which all the 
local spins are are polarized along the positive direction: 

B(^) . . . B(fi M ) \t . . . t) (4) 

Each B operator lowers the total spin of the system by one, hence the total spin of the 
state above is given by 

S = N-M (5) 

In the following, the distinguished role of the field theory ground state will be played by 
the antiferromagnetic ground state, in which N = M and the total spin is vanishing. 

The set of Bethe rapidities {/i} satisfies the system of Bethe ansatz equations, which 
for this model can be written by defining a function of the variable x: 

r>( \r -i x / \ a(x) 1 r sinh(x - fi a + ij) 2luJ 

B{x\{p), U ) = B» = ^yII sinh(x _ /XQ _, 7) ^ (6) 

relative to the eigenstate identified by the roots {//}, with 



a(x) 



7 7 M 

sinh(x — A — i— ) sinh(x + A — i— ) 

sinh(x — A + i— ) sinh(x + A + i— ) 



(7) 



d(x) 

Then the Bethe ansatz equations (BAE) are expressed like: 

B>,) = -1, j = l,...,M (8) 

It is possible from ([3]) to define a "twisted" transfer matrix, whose eigenvalue relative to 
the state |* W ({A})) is: 

rMM) = a( X ) n S ' nh( f - X > . + . i7) + e^ d{x) n ^ - - ^ (9) 
j sinn (a; — Aj) j smh(x — Aj) 

In the following, it will turn out to be useful the use of the rescaled variables (the 
"rapidities" of field theory) and inhomogeneities, and in particular of: 

9 j = % 3 (10) 
The sine-Gordon model is then recovered by the double scaling 

71 n, 4 

= -A ->■ oo , a ->■ , ae — const = — (11) 

7 m 

with m the soliton mass. 
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Because of the m-periodicity of the Bethe equation, one is free to chose to restrict the 
rapidites to the fundamental strip as — f,f , which corresponds, in terms of rescaled 

rapidities, to the strip — |(p + 1), |(p + 1) . 

The logarithmic form of ([6]), i.e., of (]§]) when considered as a function of one rapidity 
fij, defines in this region the counting function: 

Z(x) = -ilogS(-z) (12) 

7T 

which takes its name from the fact that solutions of the logarithmic form of the Bethe 
equations are characterized by: 

Z(\ J ) = 2nI J , J iGZ + ^, 6 = (N - S) mod2 (13) 

for some integer or half-integer J, with 5 = 0, 1 specifying the sector of the field theory 
the state belongs to, as explained in, e.g., [30] . 

Here we give only a short summary of the properties of the counting function, more 
can be found in [311 ESI E2J EHl ESI E3 EH US]- The function Z is real analytic and 
satisfies the nonlinear integral equation (NLIE): 

Z{6) =mL sinh0 + g{6\ {I}) - i E f dxG{6 - x a ) log FD (l + (_) V 7 ^) + a (14) 
with a = p^-uj for neutral states and tu < tt/2 and 

N H M c 

9(0\{I}) = Ex(*-Aa)-£x(*-Ao) (15) 

- E (x(« - y~) + x(« - ^ + )) - E x«(0 - A») 



X(x) = 2tt / G(x)rfx (16) 
J o 

= (p+1) /— — oos(eifc) (17) 

V ; ^ ; i 2vr cosh f sinh V ; V ; 

where we denoted x ± — x±ir), with rj some small real quantity, and the subscript 
reminds that the function x must be considered in the second determination [29j [25], [27] 
when the imaginary part of its argument exceeds min(l,p)7r, as in the case of wide roots. 
To clarify the definition of the source term g above, we need to explain that solutions of 
([13]) can be classified according to their position in the complex plane as: 

• real solutions of the equation 

1 + (-l)*e iZ(A) = (18) 

constitute the Dirac sea in the field theory limit. They will be labeled by using a 
tilde on the variable and their set by {A}. 
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• holes are real solutions of (ITS]) that are not among the Bethe roots. Their number 
isN H . 

• special solutions are real solutions of ffl8|) in which the counting function has a 
negative derivative. Their number is denoted by M5. 

• close roots are present in pairs and have imaginary part in the strip between 
— min(p, 1)tc and min(p, 1)tt. Their number is Mq- 

• wide roots are present in pairs and have the absolute value of the imaginary part 
in the strip between min(p, 1)tt and 7c(p + l)/2. There will be Mw wide roots. 

• self- conjugated roots sit on the boundary of the periodicity strip and have imaginary 
part equal to + 1). 

The ground state of the sine-Gordon model is realized as the unique state which has 
all the roots on the real axis, quantized with all consecutive half-integer quantization 
numbers. In the scaling limit, excited states are completely specified by the quantum 
numbers of the holes and of the complex roots, whose combinations correspond to the IR 
excitations. The rapidities corresponding to these objects are determined according to 
( TT3|) . and self-consistently determine the source term in f fT4|) . hence the counting function 
itself. 

The total spin of an excited state above the antiferromagnetic ground state will be 
called, in field-theoretical language, topological charge Q of the state itself. One can 
paraphrase the relation ([5]) for the total spin in terms of the excitations defined above 
as: 

Q = 2S = N H - M c - 2N S - (1 + sign(p - l))M w (19) 

In general [25], it can be stated that in the repulsive regime, each hole in the source 
terms carries a unit U(l) charge, which in the language of spin chain corresponds to a 
unit spin, i.e., to a missing creation operator. Such a charge is lowered by one for every 
close root and by two for every wide root. On the other hand, in the attractive regime, 
wide roots correspond to independent excitations, not carrying any U{\) charge. We 
then expect that the wide roots correspond to creation operators, therefore lowering the 
spin by one, only in the repulsive regime p > 1. Conversely, in the attractive regime 
p < 1, wide roots enter the expression of a state only through their effect on the other 
roots and the determination of a counting function. 

With respect to the infrared description of the spectrum, it is known [27] (see also 
[35]) that 

• soliton and antisolitons correspond to holes in the Fermi sea, quantized with half- 
integers. 

• the solitons polarization states are described by arrays of the first kind, having 
common real part and containing in any case exactly one pair of close roots. The 
position of the roots for asymptotically large volumes [29J can be determined with 
exponential precision in the size and is reported here below. Arrays can be either: 
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— odd degenerate 



O = 9 + in 

with real 9; 
even 



p + 1 



6 k = 6±iii ( — kp 



A; = 0,1,..., 



2p 



(20) 



9 k = 9±in(^-kp^ fc = 0,l,..., 



2p 



(21) 



• the breather degrees of freedom, when p > 1, are described by arrays of the second 
kind, containing wide pairs only. These are, for mL — > oo: 



- odd degenerate 

p+1 

9 = 9 + iit^— 



9k = ± in [ — — — — kp 



fc = 0,l,...,s (22) 



— even 



4 = 9 ± ivr 



k — 0, 1, . . . , s 



(23) 



with < s < ^ — 1 and real 9. In particular, they describe the 2s + 1 and the 
2s + 2 breather, respectively. 

For simplicity, we will consider in the following only states in which the number of 
special objects is null and the counting function is monotonic on the real axis, which is 
the case for sufficiently large values of the size. With proper modification, the treatment 
can be in principle extended to account also for non monotonic counting functions, but 
this appears to be more cumbersome and will not be reported here. 

The function Z{9) is suitable for numeric computation and can be determined in a 
time of the order of minutes for the simplest root configurations. For a fairly larger 
amount of time, one can determine the counting function on a suitable grid on the 
complex plane, even if only its knowledge on three contours surrounding the real axis 
and the complex roots will be needed in the following. 



3 The generating function 

A convenient method for writing connected correlation functions is by differentiation of 
a generating function. In particular, we are interested in the expectation value on the 
exact finite-size vacuum ty : 

G u (x) = (* | e-*™^*)-*™ |*o) L (24) 

in which oo is here a real number and the subscript L stands for the size of the system. 
The field in the exponent is proportional to the fraction of topological charge in the 
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interval [0, x] (as in [18]) and is realized on the lattice [16J by a string of operators acting 
on the local spin Hilbert spaces as: 

o 2m 

M*)-*(0)) = f £°f (25) 
27r j=i 

with cr 2 the usual Pauli matrix. A similar operator, the sum over projectors on the 
spin-up state, reads: 

i 2m 

Q 2m+ i = ^E(l-^) (26) 
z i=i 

A convenient representation of the exponential of f)26p was provided in [T7] in terms of 
the transfer matrices f , f u of two spin chains: one corresponding to the actual physical 
system, and the other to an analogous system in which a twist in the boundary condition 
ui had been introduced. 

We make use the results of [19J, which presents the solution of the inverse scattering 
for arbitrary inhomogeneities £j / = 1, . . . ,2N. The magnetization operator is written 
as: 

m—l / m \ — 1 

<r z m =U( A + D ) (6 " *y/2) (A - -D) (£ m - ry/2) n (^ + D ) (6 " ^/2) (27) 
z=o \Z=0 / 

which allows to write down [17] the generating function in the inhomogeneous chain as: 

m / m \ — 1 

e --Er=o-f = n^f u ((-l)'A-i7/2) n^HA-^/ 2 ) ( 28 ) 
z=o \z=i / 

We will associate here the set of rapidities {/i} to the ground state and introduce a 
complete set of eigenstates of the twisted transfer matrix on the right of this operator. 
The transfer matrices act diagonally on the respective eigenstates, so that the expression 
obtained is: 

ct ■ -t. / r -i \ \ 

{ fj T„((-l)'A- n /2|{^}) 



(*o(M) 


v — \m ~ 

e ^ 3 


*o({A*})) £ 


(*o(M) 


*o(M)) L 



with 

, rm ) = l(^({AUI^o({/i})) L | 2 , . 

U ^ (*o({A*})l*o(M)>i(*({A})|*({A})> L 1 ' 

The product on the right hand side yields, in the scaling limit, the phase: 

TT r -^ T ^((~ 1 )' A ~ !T/2j{A}aj -«(P({A}„)-P(M)) /oi \ 

l\ r ((-iyA-H/2iM) ^ e (31) 

which is reviewed in O 

What written above is then a formal decomposition of the generating function of 
connected correlation functions. Its derivatives with respect to the twist u provide a 
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form factor expansion familiar in the framework of field theory. For simplicity, we retain 
the form f[2"9"j) . in which every term acts as a generating function for the vacuum-to-state 
probabilities, and refer to it as amplitude expansion. 

The intermediate states are defined by the root structure which is encoded in the 
source term of ([ED , i.e., by the number of holes, close, ... roots and by their quantization 
numbers. The question about the completeness of the presently known solutions of the 
nonlinear integral equation is still open, to our knowledge. 

The amplitudes have the property: 

„4({7}, x, t) = e-^« A »*+^« A »- £ °)^4({7}, 0, 0) (32) 

where V({\}) is the total dressed momentum of the state \P({A}) and £({A}) its energy, 
while Sq refers to the energy of the ground state. 

We introduce here a compact notation which will be used throughout the paper: 

x^ = x ± ii] (33) 

for some small 7], denotes that the variable x is slightly shifted above or below the real 
axis. 

Explicit expressions for the exact finite volume energy and momentum of a state 
identified by a counting function Z and a given set of holes, complex and special roots 
are: 

N h N a M c M w 

V = ^2 m smn hj — ^2 ( m smn yf + m srnn u7^) ~^2 m smn c j ~ X] m s ^ nn w i 

j j j j 

"I POO 

-- d6' coshfl'S log(l + e iZ(e ' +) ) (34) 

7T J-oo 

N s Mc M w 

2 /J m cosh y s — 22 m cosh Cj — ^] m cosh Wj 

j j j 

f de' sinh e' + % log(l + e iZ{e ' +) ) (35) 

7T J-oo 

for which explicit computations can be found in [351 ESI 125]- The time shift phase comes 
from applying the double-row transfer matrix along the vertical direction along the lines 
of [291 Eg]. 

Note that the integrals in (I34|35p are well defined and finite with the convention (|33|) . 
We underline here that, since the asymptotic behaviour of the counting function for large 
values of its argument is that of an hyperbolic sine, the factor e ±lZ( - x •* goes to zero faster 
than exponentially when x — > ±oo. 

4 The result 

Suppose the ket |0)l to be the (finite volume) vacuum, which corresponds to the state in 
which all the roots lie on the real axis and are quantized by half-integers, without holes. 



£ - £buik = 111 cosh h J ~ 

j 
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The set of integers {/} defines instead an the excited state in the twisted system, as 
can be used as a starting point for solving self-consistently for the source terms in ( I14p 
and for the counting function itself. The sine-Gordon sector [33] is reproduced by the 
configuration of roots having 2S + 5 + M sc e 2Z and we will consider, for definiteness, 
half-integer quantization numbers for the rapidites, i.e. 5 = 0. Hence, in the following, 
the number of self-conjugated roots is required to be even. To have a non vanishing 
matrix element, it is moreover necessary that the total number of roots in the excited 
state is the same as that of the ground state. 
The generating function is given by 

G u (x) =N„J2 e-" (p({/}) - p °U({/}) (36) 
{/} 

where 

^=(£-4({/})j (37) 

We shall now focus on the single terms of the series A, which can be computed from 
the knowledge of the counting function of the ground state Zq and from the one for a 
generic twisted eigenstate Z\. Here below and in the following, we will denote the indexes 
relative to the holes of the excited states by h and the ones relative to the complex roots, 
generically, by a c, so that a hole solution will be denoted by A^ and a complex root by 
A c . This shouldn't generate confusion with the "close" roots, as notation will be clear 
from the context. Moreover, to shorten notations, we write the signs {c a }, with the 
convention that Choies = 1, c comp i exroo t s = — 1. We define hereby the functions: 

mL sinh^fx — y) 
^ V)= Z p {x)-Z p {y) P = A '° (38) 

A(x) = 5(x) - G{x) (39) 

and: 

/oo 
duA(x-u) log(l + e iaZo{ua) ) 
-oo 

/oo 
duA(x-u) log(l + e iaZ ^ u ) ) 

= ^^-^lQg ^^,, (40) 

With these definitions, the expression for the term associated to the twisted state 
{A} and the vacuum, evaluated in the origin, is: 

Aii}) = fl^n^n^)^ (4i) 

in which products are taken over the complex roots A c and the holes A^ that define the 
excited state and the quantities below appear. 
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The complex roots are taken into account by the factor: 



z (\ c ) 



n cos — " - — r f°° 
C(\ c ) = ry, f \\ ex P 9 / du^C + (u) coth -(A c - u) 



7 ^(A c ) 



+2/ rf-u y^c Q ,G'(A Q , — u) logsinh — (A c — u) 

J — oo „ 7T 



(42) 



where the sum runs over the holes and complex roots X a with the convention stated 
above, with the usual understanding that the second determination has to be used for 
G whenever X a is a wide root. The concerned reader may note that the logarithm of 
an hyperbolic sine grows linearly when its argument go to infinity, while the factor G 
decreases exponentially, hence the last integral is convergent. The holes enter in the 
result by the factor: 



n(\ h ) 



7T COS 



7 Z' x {X h ) 



exp 



2sign(l -p)M w - 2 / du ^ c a G(X a — u) log (f\(\h, u ) 

J-oo a 

'- 00 du 



— 3 [Cq (u)d x logip {X h ,u + ) - C^(u)d x fog(p x (X h , 



U 



(43) 



while the term V contains the various sources mixed: 



V 



Uch sinh ^(A c - X h ) sinh ^(A c - X h ) 



(44) 



The factors that embody the contribution of the "Fermi sea" of Bethe roots in the 
thermodynamic limit are denoted by $ and S, the former being: 



exp 

- E 



M 2 W 



a,a f =± 



/oo roo 
dx / dyJ2 c aCpG(Xa- x)G(Xp - y)\og(f X (x,y) - 
-co J-oo a3 

dxC°\x) [°° dy[Cl(x)dljog Vx (x°,y°') - C°(x)dl y \ogMx a ,y c 

3 J—OO 



/'OO £^<£ poo 
— G(X a -x) / dy% Ct (x) d y fog if x (x,y + ) - C£(x)d y log if (x,y^ 
-oo 7T J—oo 



again with the same conventions on the sums over sources. Note once again that the 
function fogfx is linear in x or y when these are large, therefore the kernel G ensures 
convergence of the integrals in the first line. We also have the factor: 



(4 



exp 



+ 2 E ^ / dx I dyA(x)£ a (y) 



-i ^ a ^IZ dx IZ dyC ^ x)CAy) \ 



sinh 2 ^(x a — y a ' — in) 



aa'=± 



+ 



dx 



7, 



dyA(x)A(y) logsinh — (y — x — in 



7T 



a 



7T ^ 2ni 



sinh ^(y u — x) 



sinh - (y a — x — in) sinh - (y a — x + in)- 

(46) 
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where we have made use of the notation A to denote: 



/oo A f roc 

duA(u)f(u) = c a f(X a ) ~J2c a duG(X a - u)f(u) 
-oo „, „, J —oo 



(47) 



and the sum is over all the sources, but the prime excludes the wide roots from the 
sum if p < 1. We would like to recall once again the fact that the integrations over 
variables with a superscript are not on the real axis, according to or convention ( |33l) . 
The argument of the hyperbolic cosine in ( 14TI) is given by 



E = -y2(\ a - f duuG(X a - u)] + ^1$ / — \og\ 

a V J-oo ) 2p J-oo 71 1 



X _|_ e iZx(u+) 
_|_ e iZ (u+) 



(48) 



and once again the first integration is convergent due to the asymptotic behaviour of G 
and the second due to the one of the counting function. The last term TZ is a ratio of 
determinants: 



n 



(49) 

Now we introduce the integral operators appearing in the overlap determinant. They 
depend either on two complex variables w, v and on two species indexes a, a' — ± or on 
the rapidities that define the excitations. They read: 



det 


1 - W , x 


det 


1 + (A + -1)G^'_ 


det 


1 - w> 


det 


[l + (A? - l) 


Guj 


det 


1 - W _ 


det 


1 - W x 





W( (X c ,v) 



W x ,o(X c ,X c 



G-u \X C 



ResA(X c ) 

\ _|_ e iZ (\ c ) 



^ (to, A c 



ResA(X c ) 

\ _|_ giZ (A c ) 



(£_, (A c -A c + F_ w (A c ,A c 0) 



(50) 



and 



W&{\ h ,v) = -^^(G„(x h -v°') + Ft(X h ,v)) 

w : x (w,x h ) = ± x f^ZlU ^ K - **) + K K a*)) 

Wo,A(A A ,A ft /) = - \, ^ {G u (A ft - A y ) + F M (A h , Ay)) (51) 



A(w + m/2) = A + H = e -i(^(^)-^(^))-2^ ( 52 ) 
for which another expression is provided in Section 16.11 
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The determinants are of the Fredholm type, integrals are performed on the real axis 
and the the species indexes and the excitations variables are summed over as well. In 
facts, complex roots of the bra state must be explicitly summed over as well as holes 
of the ket, if any, subtracted, as exemplified later in (192]) . The source function needs, 
for some configurations, to be evaluated in regions in which the imaginary part of the 
argument exceeds minfjo, l)ir: it is therefore necessary to use the second determination 



The integral operator G w is defined as: 



GJw) 



dk 
— < 
2vr 



: ifkw/TT 



sinh 




) k + too 




e~ iul sinh ^ 


+ sinh 




k + iu 



(53) 



and reduces to (IT7|) for u — > 0. Moreover, 



OO n QQ 

/ dx\ . . . dx n Gu (w a — Xij (1 — A(xi)j G w (x\ — x%) . . . 

n=l J -°° 

. . . (l - A(x+)) G^ (x+ - v-') 



(54) 



with and obvious extension to the case where one or both rapidities in the argument 
appear in the source. 

For what the norm determinants are concerned, their expression can be written as: 



det 



W x 



(55) 



where x stands for one of the two states and: 



cr,(T 



W x (X h ,X h/ ) 



2vr 1 + e - iaZ °( w,r ) 
1 1 

2tt 1 + 

l l 



G (w a - v a ') 
G (w a - v°') 



2vr 1 + e- iaZ ^ w ") 

ry/ ( \ \ G i^h - Aft') 



G(w a -\ h ) 



(56) 



where the "holes" terms are present for a generic excited state and 7r~ is a number 
slightly smaller than n. 

An interpretation in terms of pseudoparticles is possible: the finite-size vacuum can 
be written in terms of the fundamental excitations (solitons and antisolitons) of the 
infrared theory, which occupy the available levels according to a (complex) filling fraction 
containing the vacuum and excited pseudoenergies. Excitations constructed upon such 
a vacuum interact both among them and with the background pseudoparticles and the 
matrix elements of the operators show such features of the theory. 
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5 Scalar products and norms 



We need to perform the computation of the scalar products of the ground state with a 
generic "twisted" state [37J EH]- Let /ii, . . . ,/zm satisfy the system (jSJ), with twist oj to 
retain full generality and Aj, . . . , Am be generic complex numbers. Then 

(0| n CMMM)) = lla=1 det H({fj,}, {A}) (57) 

3- 1 n sinh(/i a - sinh(A fe - A a ) 

a>b 

in which if is 

M M 

H? k = a(\ j )t(/ji k ,\ j ) nsinh( y u z -A J -2 7 )-e- 2 ^rf(A J )t(A i ,ii fc ) JJ sinh(^-A i +2 7 ) (58) 

i=i i=i 

with 1 < j, k < M and 

—i sin 7 „ , 

t(j*, A = — — — -I — - -. 59 

smh(/i — X) smh(/i — A — ry) 

From this expression, we can extract both the overlaps and the norms of the states 
after some manipulation, in which we make explicit use that the rapidities {A} also 
satisfy Qg). 

Two alternative expressions, which are suitable for the scaling limit, can be provided 
for the overlap. We refer the reader to |T7], where the determinant of the overlap matrix 
was written as a Fredholm determinant on a contour. With some variation of their 
method, in which use of the Bethe equations and of the definition ([6]) is explicitly made 
and whose details can be found in appendix El the scalar product (1571) can be written in 
a form which is more suitable for subsequent manipulation. In the following expressions, 
the quantity uj denotes the relative twist of the state. 

e~^ M n 3 - d(\ 3 )d( N ) (1 + B^(Xj)) sinh(^ - A fc + 27) dct A fr- 
cosh(£ Xi - Y<Hi) j t l sinh(/ij - Afc) ^ 

_ Uj d(\j)d(nj) (1 + B x (nj)) sinh(A i - fi k + ij) ^ / _ ^ 



cosh(£ A z - £ Hi) sinh(/i i - A fc 



with the matrix 



jj-u K_ w (\j - Afc) Hi sinh(Aj - m) „ sinh(A j - A; - ij) 
hk 1 + 5 M (A,) n^sinh^-AO YsinKA. -^-n) 1 J 

^ = KjjMjj^JMk) Hi sinh(/Xj - A;) ^ sinh^ - ^ - nO 
1 + ^(^0 Il/^j sinh(/Xj - /i/) ^ sinh(/i J - — A { — £7) 



and the function: 



KJx) = coth-(x-z?r) - e 2iw coth -(x + ivr) , AT (x) = AT(x) (63) 

7T 7T 
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For the computation of norms, one considers the limit {A}, {p} — > {y}, for which the 
matrix above becomes simply: 



^1 = ? { Vf k \ (64) 

For the remaining part of the section, we shall be using rescaled rapidity variables. 
Moreover, unless otherwise specified, we shall consider the state \i to be the ground state 
of the (untwisted) inhomogeneous chain, while the state {A} is considered as having a 
twist. 

As a preliminary step, one observes that by applying the definition of counting func- 
tion and a representation of the cosine as an infinite product: 

i+w=^^> nji-^y (66) 

from which, considering the ground state with 2M roots {//}, having labels ranging on 
half-integers between — M + 1/2 and M — 1/2, we have: 

A) JL/A z M) ;r ( M + i-^)r(M + i + ^) 

(66) 

The last ratio tends to unity in the limit in which the number of roots goes to infinity 
and will not be rewritten in the following. The case in which there is a finite number of 
excitations yields the same result, if one considers in the product above the set {//} of 
all the real roots of (fl8|) . 

Having established this fact, we consider the state {A} to be excited and the state 
{p} to be the vacuum, identified by the subscript 0. We are moreover interested in the 
normalized matrix elements, so we divide the overlap by the norm of the two states. 

Let us multiply and divide by the holes and the complex roots, in order to obtain 
expressions in which all and only the real solutions appear. This is convenient in that 
we can consider the ratio between each hyperbolic sine appearing in the denominator of 
the expressions ( |60|) and the differences of the counting function computed at the points 
in the argument of the sine, as arising from the product representation f l66l) . Following 
[TT] this defines the functions ( 1381) . 

After illustrating the general procedure, it is simpler to consider two additional states, 
whose rapidities we label by {p}, {z/}: at the end of the computations, we will send 
{p} — > {A} and {z/} — > {p} and show that the poles arising from the factors of the kind 
( I18p are canceled by the zeros of the hyperbolic sines in the expression for the scalar 
product. In order to obtain a product involving only the real solutions, we consider the 
ratio: 

sinh 3-(Aj - p k ) sinh ^(fij - v k ) = IL C ' sinh ^(A e - p c >) Uh,h' sinh ^(A fe - p h ,) 
j,l sinh l^i ~ Vk ) sinh ~ P k ^ lUsinh l(\ c - p h ) sinh ^(p c - \ h ) 

j-j- sinh - A fc ) sinh - v k ) sinh l(\ h - p,j) sinh ^(p h - p,j) 

j,k sinh - (j, k ) sinh Zfo - p k ) h j sinh ^(A fc - pj) sinh l(p h - ~Xj) 
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sinh l(\ c - pj) sinh ^(p c - A 3 -) 

On the other hand, for what the factor involving the counting function is concerned, we 
can write: 

1 + e iMpj) i _|_ e ^o(Pc) x + e ^o(p h ) /! + e iz A ( Ph ) \ 2 i + e iZofe) 

AA 1 4. e iZx(j>j) ~ AA 1 _|_ e iZ\{p c ) AA 1 4- e iZ\{p h ) AA I i _|_ e i^o(Pfe) J AA ]_ _|_ e iZ x {pj) ^ > 
3 c h h \ '3 

By multiplying the first and the second product with the first ratio of ( 1671) and taking 
the limit to coinciding states, we obtain the term (144jl . which already contains a finite 
number of rapidities, apart from a phase factor. 

The factors in the first term that contain the same index for the hole, together with 
the third term of (|67|) and the third term of (168]) yield: 



l\U{\ h ) , W(A(/i)) = 11 — 7T rr^ (69) 

while the last of (1681) . accompanied by the corresponding product in the rapidities, 
with the second in (1671) provide the factor: 

$ _ V ? o(^j^fc)v?A(Aj, Afc) ^ 7Q ^ 
Vo(a*j> Afc)^A(A*i, A fc ) 

From all the previous expressions, one also obtains a phase factor containing sum over 
rapidities of the difference of the two counting functions, that will be of no relevance in 
the following. 

The last product in (l67j) is already in a form suitable for the scaling limit; together 
with the part of the product in the first term that contains the same index for the close 
root, it may be rewritten as: 

and constitutes a multiplicative contribution from complex roots. 
According to our previous analysis, we write: 




l(*({A})|*(M)>| 2 



,. , ( , -sv<t> l[U(\ h )l[C(\ ( 



ii*(M)ini*({A})ii 

Where the definitions (1441) and 



det (l - U^) det (l - U, 



det (1 - U x ) det (1 - Urf 



(72) 



s = tt sinh(/jj - A fc - ij) sinh(Aj - p k - 27) 
j£ sinh(/i i - ji k - ij) sinh(A i - X k -ij) 



have been used. 
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6 The scaling limit 



Assuming that a function / is analytic in a connected region of the complex plane, that 
contains all the roots {A}, satisfying (jSJ), then the sum of the values of the function when 
evaluated in this points can be written as 

^ /(Aj) = Loo 2^ (cosh( M - 6) + cosh( M + 6)) 

N h M c M w M s 

- £(A* f)(h k ) + £(A * f)(c k ) + ^2(A n * f)(w k ) + 2 ]T(A * /)( Sfc ) 
fc=i fc=i fc=i fc=i 

+*5> /_^(A^)K)log^(l + e^)) (74) 

where the symbol * represents convolution over the real axis and the distribution (|39|) 
appears. By using this expression, which can be obtained in an analogue way to the 
NLIE itself [321 1251 1281 127] , one can compute the scaling limit of the factors above. 

6.1 The determinants 

We want to reduce the expressions above to a standard Fredholm determinant form 
det(l + K) = e^™^ -1 )™ ltT l Kn ]/ n , p or an analytic function (at least on the real axis), one 
can write 

y /(A,) _ r dw f(w) 

yih.-M:, A;) J2mUi(w-X l ) 1 ; 

This is applied to the matrices f )106p . by considering the variables A Q , A^ as two complex 
variables w, v integrated on a closed contour. In order to do so, one employs the function: 

. . il sinh-(w — rii) sinh ^ (w — A/ — m) , , 

A (w) = II . hX( n — ; \, 76 

sinn ±(w — Hi — iix) sinn - L [w — Ai) 

The latter expression contains both the poles in the values of the A roots to be summed 
over and the zeros in the values of the \i roots to be avoided, a fact that allows us to 
keep the contour of integration under control for every couple of states. Be the reader 
aware that we are using rescaled variables. 

If we consider first the state {//} to be associated with a (twisted) excited state and 
the state {A} with the ground state, it is sufficient to consider a contour that encircles 
the real axis. We underline that all the zeros of the factor 

l + 5^(A a ) -> 1 + (77) 

are all the real roots and holes in the state {/i}. This means that the product 

sinh-(w - 

i * 

will cancel all the poles corresponding to real roots, but not the ones corresponding to 
holes, which will be treated separately. 
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With w a generic complex variable, having < \^sw\ < 7rmin(l,p) strictly, it is 
possible to exponentiate the product and apply the formula (173j) to the sum of logarithms. 
It is however necessary to choose the contour in a way to avoid the branch cuts, which is 
simply done, also numerically, provided 7 is not too close to tt or to zero. The result is: 



A(w) 



n c sinh £ (to - A c - iir) l\ h sinh £| 



w 



Uc sinh I (w - A c ) lift sinh £ (tu - A h 
exp 



itt) 



TT 



cr=± 



7, 



itt) — coth — (w 

TT 



+ / dx c a G(A a — x) log 



sinh 



X — W — Z7T 



sinh 



2, 



x — w) 



C a (x) 
(78) 



which holds whenever w is not on the real axis. However, it is most simply written if one 
shifts the argument by in/2, which yields ( 1521) and can be used to compute the factors 
needed in (15T1) . provided p ^ 1/2. The form above, instead, seems to be more useful for 
numerics and for extracting the expression ResA(X c ) in (150]) . 

If the counting function in ( J77|) refers to an excited state, then it is necessary to 
subtract from the sum over poles the unwanted ones corresponding to holes, a task 
which is performed by using a term like: 



ly 

tt , , sinh - 

holes 7r 



(w - A, 



'ft, 



v) 



(79) 



as long as the number of holes is of order of unity. 

It is also possible to extend our analysis to the case of {fi} being the ground state 
and {A} an excited state: holes need not to be subtracted anymore (the ground state 
has all the Dirac sea filled), but complex roots outside the contour need to be explicitly 
added when computing the Fredholm determinant. 

This cannot always be done by deforming contours, because of the poles in the kernel 
at it — v — ±7r, ±irp. It then follows that roots lying beyond min( : |, 22) must be treated 
separately and enclosed in different contours. 

Taking into consideration neutral states with rapidities quantized with half-integers 
[2"7j S3] the close pairs of complex roots approach their infinite-volume position ( 1211) 
keeping their distance larger than tt. Moreover, as was observed in Section [2j the wide 
roots do not correspond to any creation operator in the attractive regime, since their 
presence does not modify the total spin. Then, in the attractive regime, the prescription 
for the external contours is to surround the region whose imaginary part is | < 1 5s9 \ < Trp. 
On the other hand, it is known that the antisymmetric soliton-antisoliton states are 
described by a pair of close roots, so once again, the prescription applies. 

The Fredholm determinants are computed on a contour: 



det 



i-lo 

TT 



6XP I ? \ 6)7 2^ • • • / ^0(vi,v 2 )6(v 2 ,v 3 ) . . . 6(v n , Vl ) 



(80) 

The contours have to surround all the roots on the real axis; moreover, for excited states, 
they also need to encircle the complex roots. In principle, one can surround each root by 
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Figure 1: contours surrounding roots of the Fermi sea on the real axis and complex roots 
for p < 1. The dotted lines are at = ±|, ±np. 



a small circle, taking care to avoid that two points in contours are separated by 
This is the prescription for the repulsive regime. 

In the attractive regime, the contours encircle only close roots, with \Q6\ < n p. In 
particular, it is known that the pair of close roots describing the polarization of a soliton- 
antisoliton pair, has an imaginary part which reaches the values =ki~ from above. The 
same arguments extend to all the close roots quantized with half-integers. Hence, our 
contours surround the region of the complex plane ~ < \Q8\ < np. 

Indicating by a subscript the states that enter in the matrix element, the overlap 
kernels, as function of the field theory (rescaled) rapidities, are: 

U ^ V) = 2ri TT^ Ku{x ~ V) '^ sinh(^(x - X h )) 2nZ>(\ h ) K " iXh ~ v) 

(81) 

and the norm kernels 

1 1 

UJx.y) = .„ , , K(x — y) 

U ^ V) = ^TT^ K{X ~ V) ~l^ sinh(^ - X h )) 2, Z \{X h ) K{Xh - V) 

(82) 

The previous expressions hold for w not on the real axis. They are valid also in the 
scaling limit but, if the limit is taken at this stage, the integrations with the kernels 
above are diverging. Therefore, we need to manipulate further these expressions, which 
will be done in the following section. 
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6.2 More on determinants 



The integral operators in (I81|l82p above are not of trace class when the integration is 
performed over the whole real axis. To proceed with the numerical evaluation, one 
would be forced to introduce a cutoff; the function A above, which is present in the 
overlap integrals, tends exponentially to unit value when its argument has large real 
part, for any couple of the states. It follows that the asymptotic behaviour for large 
rapidities of the "overlap " kernels is the same. Then, by considering normalized overlaps 
as in fl30|) . one may argue that the result is independent from the cutoff and the ratio of 
determinant such as flHUj) is finite. 

However, we will not need to follow this path. Instead, we start from integrations 
which surround the real roots of an eigenstate of the finite chain, so that contour integra- 
tions must surround only a finite number of roots. Then, we split integrations and group 
them into series. Finally, the boundaries of integration are sent to infinity, providing an 
exact summation of these families of terms when the scaling limit is taken. As a result, 
we retain a Fredholm determinant form, but with a trace class kernel. 

Considering the kernel in the expression of the norm, the integral on the contour can 
be written as: 

dw K{w — v ) f A dw ( K{w + — v) i K(w~ — v) + .A 



2ni 1 + e iZ ^ J-A 2m \ 1 + e ~ iZ ^ w+ 1 1 + e iZ ^ w ) 



where the integration boundaries are intended to be more negative and more positive 
than the real roots with largest absolute values, and the parts of the contour which are 
perpendicular to the real axis are assumed to give negligible contribution. By writing 
down explicitly the first terms of the Fredholm determinant expansion with this form for 
the kernel and noting that the function K has no poles on the real axis, it is possible to 
show (see appendix IB]) that: 



V- </"—!- / dw n K(wi - w 2 ) K(w n - wi) ^ 



?? 



%K 
~ 2^ 



1 Jdwi 7dw n ( G{w?-w?) G{w?-w?) \ 

A. / 2tt '"i 2tt { l + e -^ CT1 ) "' 1 + e^ z ^) ) {Q > 



T) 

n a 1 ...a n =±_ 



for some generic counting function Z, where we have used the fact that the kernel in the 
NLIE ( Tl4l) can be expressed [16] in terms of the spin chain kernel as: 

00 °r -ik(a\ 

G{w) = 52(-lT-\-iKr(w) = / dqe™* .™ (85) 

n=i _4 1 - iK(q) 

in which the superscript *n denotes n-times convolution. 

The first term in the right-hand side of ( }84"1) is formally divergent. However, it can 
be combined with the analogous terms arising from the overlap determinant and written 
as a series, as seen in Section 16.41 
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In order to deal with the overlap determinant, one can similarly set: 

-iKM 



n-l. 



-tK w y n (w) 



n=l 



dqe 



iwq 



(86) 



The Fourier transform is performed by separating the asymptotic behaviour of the func- 
tion and by regularizing through a damping exponential. A reasoning analogous to that 
for the norm leads to the following overlap determinant: 



det 



1- E W^{w,v) 

cr,a'=± 



(87) 



where we have defined a function W of two real arguments and with two contour sub- 
scripts. To be clear, we denote, e.g., W+_(x, y) the expression resulting from the dressing 
of the kernel of the overlap determinant, where the first variable has real part x and is 
on the upper branch of the contour and the second has real part y and is on the lower 
one. Explicitly: 



A{w a ) 



27r V 1 + e~ aiZ ^'' 



^-6 a , + (A (w + )-lfjG^(vf-v) 



(88) 



which is manifestly of trace class due to the presence of the soliton and antisoliton filling 
fractions and to the asymptotic behaviour of the function A, which tends to one for large 
values of the argument. 

We will shift the contours of integration to a definite value = rj and define the 
kernel W with real argument only. 



W%{w,v) 
W( (X c ,v) 
W^(w,X c ) 



1 / A(w + iarj) 

27T I 1 + e -icZ {w+i<jri) 

ResA(\ c ) 

\ _|_ giZo(A c ) 

1 / A(w + iarj) 



S a+ {A + {w) - l))G- u (w-v + i(a - a')r,) 



2n \ 1 + e~ iaZ °^ w+iar ^ 
( ResA(X c ) 



X _|_ e iZ (\ 



G-u (X c -v- ia'r)) 

- 5 a+ (A + (w) - l)j G-u (w - A c + iarj) 

1 J G-u (A c - A c /) (89) 



W ° x (w,\ h ) 
W , x (\ h ,\ h ,) 



— tTTTT w (^h-v- ia rj) 

I ( i^ra " ^ te» - 1)) <M«, - « + - 

(A-\\ h " ^ 



-l)G u (\ h - A,0 



(90) 
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It is then possible to apply again the partial summation of some families of terms in the 
Fredholm determinant series along the lines of appendix [Bj In this case, the terms that 
can be summed over are the ones with domain only in the upper branch, i.e., those like 



2tt 



-<U AT 1 



w 



i a 



(91) 



Unfortunately, it is not as easy as in the previous case to Fourier transform and sum 
over all the terms, hence the result is given in the rather involved form f l54l ) . 

The expressions above imply explicit subtraction of holes in the ket and summation 
of the complex roots of the bra, when used to compute the determinants of ( 14"9"|) . For 
instance, when multiplying W\$ with itself, one has: 

/ \ H /"oo 

(W 2 xfi ' (w,v)= £ / dzW^(w,z)W^-(z,v) + J2W+ (w,X c )W^(X c ,v) (92) 



and similarly for Wo t \. Then one can make use of 



det 



1 - W 



OO 1 



ra=l 



n 



W r 



(93) 



with the above recipe for the product and the trace. 

At this stage, contour indexes can be interpreted as species indexes, as they enter in a 
symmetric way in the result. In order to do so, one could deform the contour up to \Qw = 
ir/2\ and, following [32], define the (complex) soliton and antisoliton pseudoenergies as: 



-(d) 
(d) 



-iZ(6 + i 



71 , 



iZ{6 



.7I\ 



(94) 



Then we understand the factor containing the exponentiated counting function as a 
filling fraction. This interpretation is very suggestive of a possible extension to two- 
point functions of the validity of some sort of Leclair-Mussardo formula [8], provided 
suitable excited state pseudoenergies are used. To avoid heavier notation, we write the 
final result in terms of the counting function alone. 



6.3 The prefactor 

As a warm-up, we look at the argument of the cosine in (1601) . By applying the formula 
flT4")) . we see that the term proportional to N cancel in the difference and we obtain 

(Xj - Hj) = - c a ( X a - ( duuG(X a - u)j 

j a \ J-oo J 

_y / a _ / d G f a _ w ) l g' iX^^ (95) 

^ ± Loo2m\ J V ' J l + e ^(« CT ) V ' 

with X a being, as usual, the roots that define the source terms in ( fl4l) . Then, we 
integrate by parts the second term and arrive to the expression ( )48l) . using the asymptotic 
behaviour x(±oo) = ±| ^ [29]. 
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We now turn to the analysis of the factor (1731) . Here the complication lies in the 
double product, but there are no conceptual difficulties in exponentiating this expression 
and performing the sum of the resulting logarithms with the procedure described above, 
since the arguments of the logarithms never cross the cut. The first step is 

L sinh(Aj - /i fc - ij) ^ sinh(/Xj - X k - 
• sinh(Aj — Afc — ij) sinh(/Xj — /ifc — ij) J 

sinh(Aj — x a — £7) 
sinh(//j — x G — ij) 




j 

00 



I duA(u)log Sm ^~ " (96) 
J -00 smh(yUj — u — 27) 

where both counting functions are relative to a finite number of rapidities. Then one is 
to apply again (JTijl to the j index. After passing to the scaling limit (TIT]) we find the 
result (146]) . with the counting functions satisfying (1T4"]) . 

The complex roots are, by definition, away from the real axis. Since in the product 
defining the complex root factor (|7TT) there appear all and only the real solutions, it is 
natural to consider the logarithm and chose a suitable contour around the real axis to 
perform the sum over roots. The resulting expression is given in (142 p . 

We turn to the analysis of the factor containing the position of the holes. The 
counting function itself may be non monotonic in some region of the real axis for some 
class of states and at small volumes mL ~ 0(1), a circumstance which is connected with 
the appearance of special roots when the counting function also crosses a quantization 
point within that region. We assume that this is not the case, even if, in principle, it 
can be worked around by splitting the sum over different regions in which the counting 
function is monotonically increasing. Having written the product (169]) in terms of positive 
functions, we can now take the logarithm and perform the scaling limit. 

Such limit is somewhat simplified by the fact that the sum is performed only over 
real solutions of the (ITS]) . It can be performed by standard techniques, but some details 
are in order. After applying the summation procedure, one of the terms has the form 

x 0= f° ( N + N ) i og (97) 

J-00 2ir \cosh(u — Q) cosh(w + 6) / (f\(x,u) 

This integral is similar to the ones needed for the computation of energy and momentum 
([25], see also [35]): one uses the fact that the counting function of a finite chain has a 
well defined limit when the argument is sent to infinity 




A% + gvrS + 27rsign(p - l)M Wi + 2u , , 

-A% - 2^tt S - 27rsign(p - l)M wt + 2u ' ' ^ ' 



where M w ^ and M w ^ are for the number of wide roots below and above the real axis. It 
follows that: 

X = signQo - l)M w (99) 
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This integral enters both in the evaluation of the "hole" factor and in the product ( 1701 . 
Here below, we report the result of the summation on the first index: 



y j log — -pz ^ — = 

j,k V 9 A(Ai,/ifc)v?o(Aj,/iA : ) 

du / l + e i^Z x (un roo ¥o(fij,X u 



( roo j_ _|_ v." I roo 

+ J2c a duG(X a - u) (log (p x (Xj : u~) - log^o(A } ( 10 °) 

the second step, again by taking into account the asymptotic behaviour ( |98|) . can be 
analogously performed and yields (145]) result. 

6.4 A note on normalization 

The "normalized" amplitudes A are defined from the bare amplitudes ( 1301) by dividing 
by the factor 

= det[l + ^]det[l + ^] 

det[l + K] 2 V ; 

with kernel (|63l) and a; on the real axis. The operators are not of trace class and we 
have not found a way of directly evaluating these determinants up to present. Moreover, 
at first sight, it appears that the factor fllOip can be written in terms of a difference of 
Fredholm series, whose terms are one by one divergent. 

However, a definite expression for this factor, in the form of a series, can be given 
considering (156]) for x = 0, for which the operator f[2"4"]) reduces to the expectation value 
of the identity on the finite-size vacuum. Therefore, imposing 

= 1 (102) 



for every size L, we obtain the exact normalization in the form ([3] 

Let us add that, being the series for the generating function derived from 
expect similar convergence properties for (1361) directly. However, it seems that this issue 
should be tackled by extensive numerical analysis. 



7 Conclusions 

We have presented an exact expression for the generating function of connected corre- 
lation functions on a cylinder, where the compactified direction is space, for the sine- 
Gordon quantum field theory. 

To take into account corrections which enter as exponentials in the size, the knowledge 
of the spectrum and of the form factors of the theory in infinite volume is not enough. 
To circumvent this problem, the computations were carried on in the framework of the 
Destri-De Vega lattice regularization, formulated in terms of an inhomogeneous XXZ 
spin chain. 
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In this framework, the problem was similar to the study of the expectation value of 
the magnetization of the spin chain in a given interval: due to the available results for 
this model, the most important of which is the solution of the quantum inverse scattering, 
we managed to write the exact vacuum expectation value of the generating function in 
the form of a series in ( 136]) and performed the appropriate scaling limit of each term to 
obtain the result (|4T|) . 

Each term in our expansion is associated to one of the states of the field theory 
in finite volume and can be interpreted as a generating function of form factors. The 
determination of the exact states relies on the ability of solving self-consistently the 
Destri-De Vega nonlinear integral equation for all the allowed source terms, which is 
in general a difficult task and limits for the moment the practical applicability of the 
method to some classes of finite volume states. 

Having shown its relevance in the computation of correlation functions, we hope to 
be able to extend the analysis of the nonlinear integral equation in the future. Further 
work would be moreover required to explore the advantages and limitations of this kind 
of formalism in the actual computation of correlation functions, as well as to identify the 
features which may be generally valid for other integrable field theories. 
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A Manipulation of the matrix elements 

To compute the determinant, we multiply the matrix if by a conveniently defined matrix 
M and its inverse. In the case of a twisted state defined by the roots {A} and the state 
defined by the roots {/i}, we consider the matrix 

coshO-,- - A fc ) n^fc sinhOj - A*) 
UijLj smlif^ - m) 

whose elements are i7r-antiperiodic functions of the rapidity /ij. Its determinant is 

Then the matrix product 

Hj n M n k can be computed by considering the integral 

f dw sinh(-i7) coshQ - A fc ) n^fcsinh(w - A;) 

J 2ni smh(w — Xj) sinh(w — Xj ± ij) 17/ sinh(u> — /!/) 
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which is vanishing when the contour of integration surrounds the real axis and the strip 
[— 7r/2,7r/2] along the imaginary axis. 

Then the result of the matrix multiplication is 

[H" ■ M\ . = (-e-^)d(Xa) (1 + B«M) Y[ sinh(A, - ^ + ^ n^sinh^ - ^) 



Hi sinh(/i a - \i 



Hi sinh(/i a - A;) 
sinh(/i a - hi - 17) 



1 + B%((i a ) Ui^a sinh(/i a - hi) sinh(/i a - A/ - 17) 

(coth(// a - /if, - 17) - e 2lw coth(// a - /i 6 + 27))} (106) 



from which the Fredholm determinant in the limit in which the size of the matrices goes 
to infinity can be recovered. 



B Dressing of the kernels 

We start from a situation in which the number of roots is finite. Given an integral 
operator Q and the kernel K, we have that: 



det [1 - (Q - K)\ = det 



1-Q* 



l + K 



det [1 + K\ 



(107) 



Here, the symbol * is used to denote a convolution in which the extrema of integration are 
not to infinity yet they are large enough to contain all the real roots. The equality can 
be checked by taking the logarithm of the above expression, expanding and reordering 
terms: 

trflog [1 - (Q - K)\ ] = tr\K - -K* 2 + -K* 3 - -K* 4 + -K* 5 - ... 

2 3 4 5 

- [Q - Q * K + Q * K* 2 - Q * K* 3 + Q * K* 4 . . .) 
— (Q* 2 - 2Q* 2 *K + 2Q* 2 * K* 2 + (Q* K)* 2 - 2(Q * K)* 2 * K . . .) 
_i (^g*3 _ 3 q*3 ^ K + 3 g*3 *k* 2 + 3Q*K*Q* 2 *K + ...) 

_I (Q* 4 -4Q* 4 *K + ... 



= tr 



log (l - Q * (l - K + K* 2 - K 



*3 



tr [log(l + AT)] 



(108) 



which is used both for the norm and for the overlap kernel, with different Qs and either 
K w or Kq. Obviously, the symbol K* 3 means K * K * K and so on. 

After organizing the summation into series, it is safe to take the scaling limit in 
the first term of the last expression, since it is of trace class on the real axis. Then 
the * becomes a convolution over the whole real axis (*). For what the second term is 
concerned, this becomes formally divergent as the integration boundaries go to infinity. 
Nevertheless, these terms are independent of the state and can be factorized out of the 
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sum ( 1361) . grouped together and represented as a series, as explained in Section EU In 
other words, there is no need to evaluate them. 

To perform the Fourier transform of K u , it is simpler to use the Bethe rapidities and 
write the kernel as: 

K U (X) = e iuJ [cosuK Q (X) -tsmuu(K + (X) + 2tanh(A))] (109) 

where 

K+(\) = coth(A - 17) + coth(A + ij)-2 tanh(A) (110) 
By applying the residue theorem, one can see that 

sinhf^ — 7) A cosh ( ~ — 7) k — 1 

KJk) = % V2 J J , K+(k) = -i ^ 4 (111) 

V ; sinh^ +V ; sinh^ V ; 

Moreover, we write: 

dXe~ iXk 2 tanh A = -2 lim / d\(e iXk - e~ lAfc ) tanh Xe~ ak (112) 

-00 Q^OJo 

then we make use of the integral representation 

valid for > 0, to write the integral above as a sum of digamma functions. After 
further massaging and the limit a — > 0, we arrive at the expression 

cosh ( ( ~ — 7) k + ioj) 
KM = ie«" U2 . - ! k J - (114) 



sinh^ 



from which (I53p follows. 



C Translations 

Within the formalism of the light-cone lattice (see [THl [36J), the operators generating 
translations along the space and time directions can be constructed by successive appli- 
cation of the inhomogeneous transfer matrix at special values of the argument: 

e***-* |¥({A})> = U L |*({A})> = f(0) |*({A})> = r(0|{A}) |*({A})> (115) 

e -^+P) mX })) = U R |M>({A})> = f(-0)t |*({A})> = r(-0|{A})* |¥({A})> (116) 

Considering two neighboring sites with opposite values of the real part of the inhomo- 
geneity: 

e-^(0|{A})e— r w (-0|{A}) ^ a(n{ x M -vm)) 



r o (0|{/i}) r (e|M) 
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e -ia(V(i\M-nW)) (117) 



By looking then at the full product in ( |29|) . one recovers the expression of momenta 
analyzed in [35J, in presence of a twist and reported in f[3~4|) . An analogous path can be 
followed by applying the transfer matrices 

f(0)f(-0) t (118) 

and performing the scaling limit on their eigenvalues, in order to obtain translations 
along the time direction. 
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